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ABSTRACT 

We address the issue of the cosmological bias between matter and galaxy distributions, 
looking at dark-matter haloes as a first step to characterize galaxy clustering. Starting 
from the linear density field at high redshift, we follow the centre of mass trajectory of 
the material that will form each halo at late times (proto-halo). We adopt a fluid- like 
description for the evolution of perturbations in the proto-halo distribution, which is 
coupled to the matter density field via gravity. We present analytical solutions for the 
density and velocity fields, in the context of renormalized perturbation theory. We 
start from the linear solution, then compute one-loop corrections for the propagator 
and the power spectrum. Finally we analytically resum the propagator and we use a 
suitable extension of the time-renormalization-group method (jPietro ni 2008) to resum 
the power spectrum. For halo masses M < 1O 14 /i -1 M0 our results at z = are in 
good agreement with N-body simulations. Our model is able to predict the halo-matter 
cross spectrum with an accuracy of 5 per cent up to k « 0.1 h Mpc -1 approaching 
the requirements of future galaxy redshift surveys. 

Key words: cosmology: theory, large-scale structure of Universe- galaxies: haloes - 
methods: analytical, N-body simulations. 



1 INTRODUCTION 

Redshift surveys have shown that the clustering properties 
of galaxies strongly depend on their l uminosity, color and 
morphological (or spectral) type (e.g. iNorberg et alj |2002| ; 
IZehavi et ai1l2004l h This indicates that galaxies do not per- 
fectly trace the distribution of the underlying dark matter, 
a phenomenon commonly referred to as 'galaxy biasing'. Its 
origin lies in the details of the galaxy formation process 
which is shaped by the interplay between complex hydro- 
dynamical and radiative processes together with the dark- 
matter driven formation of the large-scale structure. 

Attempts to infer cosmological parameters from galaxy 
clustering studies are severely hampered by galaxy biasing. 
A number of theoretical arguments and the outcome of nu- 
merical simulations both suggest that, on sufficiently large 
scales, the power spectra of galaxies and matter should be 
proportional to each other: P g = h\ P m where the linear bias 
factor &i depen ds on galaxy type but is gene rally scale in- 
dependent (e.g. Cole si 1 19931 : Im arm et al.lll99ct) . Similarly, to 
model higher-order statistics, such as the galaxy bispectrum, 
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it is generally assumed that galaxy biasing is a local process 



such that 5 



bi5 m + &2f> m /2 + . . . where 8 g and 5 m are 



the (smooth ed) galaxy and dark-m atter density contrast, 
respectively (|Frv fc Gazta naga 1993i). However, the reliance 
of these phenomenological approximations limits cosmolog- 
ical studies to very large scales whereas data with better 
signal-to-noise ratio are already available on much smaller 
scales. Mor eover, future studies of baryonic acous tic oscil- 
lations (e.g. ISugivama|[l995l : lEisenstein fe Hulll998l l will re- 
quire measurements of the matter power-spectrum with per- 
cent or even sub-percent accuracy in order to shed new light 
on the source of cosmic acceleration. Understanding and 
controlling the effects of galaxy biasing with this precision 
will be challenging. All this provides a very strong motiva- 
tion for developing more accurate (and physically driven) 
models of galaxy biasing. 

A number of authors have used the power spectrum 
statistics to explore the scale depende nce of galaxy bi- 
asing based on numerical simulations llCole et al.l 2005 ; 



Seo fc Eisensteinll2005l:lHuff et al.ll2007l: iManera et al.1 



2010; 



Manera fc Gaztafiag 3 120091: iMontesano et al.1 l2010ll or an- 



alytical calculations (iSeliakl 120011: ISchulz fc White! 120061 : 
iGuzik et ai]|2007l : [Smith et al.ll2007h stemming from either 
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perturbation theory or the halo mo del for the large-scale 
structure fsee lCoorav fc Shethll2002l for a review). The gen- 
eral picture is that galaxy biasing is expected to be scale de- 
pendent (i.e. -Pg(fc) = b(k) 2 Pm(k)) and the functional form 
of fe(fc) can sensibly depend on the selected tracer of the 
large-scale structure. 

Since galaxies are expected to form within dark-matter 
haloes, understanding the clustering properties of the haloes 
is a key step to accurately model galaxy biasing. This 
is a much simpler problem, considering that dark-matter 
haloes form under the sole action of gravity. It is in fact 
expected that long-wavelength density fluctuations modu- 
late halo formation by modifying the c ollapse time of local- 
ized short-waveleng th density peaks jBardeen et al.l 1 19861 ; 
ICole fc Kaiseilli989h . This argument (known as the peak- 
background split) predicts that, on large scales, the halo 
overdensity 5h = b £ m where the bia s coefficient b varies 
with the halo mass l|Mo fc White! 1 9961 The numerical value 
of the bias coefficient is determined by two different occur- 
rences: first, haloes form out of highly biased regions in the 
linear density field (Kaiser 1984; Porciani et al.lTl999h and, 
second, they move over time as they are accelerated towards 
the densest regions of the large-scale structure by gravity 
|Mo fc Whitelll996l ). These two phenomena generally go un- 
der the name of "Lagrangian bi asing" and "Lagrang ian to 
Eulerian passage", respectively. iMo fc White! ll 19961 ) dealt 
with the second problem by assuming that long-wavelength 
density perturbations evolve according to the spherical top- 
hat model. A more sophisticated gener alization of the peak - 
background split has been presented bv lCatelan et alj (ll998T ) 
who assumed that also the large-scale motion of the density 
"peaks" is fully determined by the long-wavelength com- 
ponent of the density field. Since the halo population and 
the matter feel the same large-scale gravitational potential, 
their density fluctuations are strongly coupled and their time 
evolution must be solved simultaneously. This makes the 
process of halo biasing non-linear and non-local even if one 
starts from a linear and local L agrangian biasing scheme 
l|Catelan et al.ll 19981 ; lMatsubarall2008l ). The bispectrum can 
be used to test this model against the standard Eulerian 
local biasing scheme (|Catelan et al.l [2000t ). 

In this paper, we present a novel and very promising 
approach to model the cl ustering of dark mat ter haloes. 
Adopting the formalism bv lCatelan et all |l998) combined 
with a non-local Lagrangian biasing scheme for the haloes 
l|Matsubaralll999h . we simultaneously follow the growth of 
perturbations in the matter and in the halo distribution over 
cosmic time. We present perturbative solutions for the cor- 
responding overdensity and velocity fields and we are able to 
resum the perturbative series in the limit of large wavenum- 
bers. Moreover, we write down a system of equations for the 
power spectra P m and Ph using th e time -renormalization- 
group (TRG) approach bv lPietronil (|2008l ) and numerically 
integrate them. Our results are in excellent agreement with 
the output of a high-resolution N-body simulation, showing 
an improvement over linear theory, and we are able to pre- 
dict the matter-halo cross spectrum with a precision within 
5 per cent for fc < 0.15 ft Mpc -1 . 

Related work has been very recently presented by 
iDesiacques et al.l (|2010f ) who computed the two-point cor- 
relation function of linear density peaks and followed its 
time evolution assuming that peaks move according to the 



Zel'dovich approximation. For massive haloes this results in 
a scale-dependent bias (with variations of ~ 5 per cent) on 
the scales relevant for baryonic-oscillation studies. Contrary 
to their approach, we do not deal with a point process but 
describe large-scale fluctuations in the distribution of dark- 
matter haloes as perturbations in a continuous fluid. On the 
other hand, we account for the full gravitational motion of 
the haloes and do not rely on simplified dynamical models 
as like as the Zel'dovich approximation. 

The structure of our paper is as follows. In Section [2] 
we present our model for the joint evolution of the mat- 
ter and halo power spectra. The initial conditions for our 
evolutionary equations are discussed in Section [3] The so- 
lution of the linearized equations is presented in Section [4] 
where we also quantify the importance of the halo veloc- 
ity bias. Using a perturbative technique, in Section [5] we 
compute analytic solutions for the propagator of perturba- 
tions (the two-time correlator). We derive 1-loop corrections 
and, in the limit of large wavenumbers, the fully resummed 
propagator. The discussion in Sections 15.11 and 15.21 is very 
technical and the less experienced readers can safely skip it 
without compromising understanding of the remainder. In 
Section|6] we numerically integrate the full equations for the 
evolution of halo and matter power spectra in the TRG for- 
malism. We then compare the results against the outcome 
of a high-resolution N-body simulation. Finally, in Section 
[7] we conclude. 



2 THE MODEL 

2.1 Dynamics of gravitational instability 

The large-scale structure observed today in the universe is 
believed to be the result of gravitational amplification of pri- 
mordial fluctuations, caused by the interaction among cold 
dark matter (CDM) particles. If we denote 8 m as the matter 
density contrast and v as the velocity, the Eulerian dynamics 
of a system of such particles, which interact only via gravity, 
is ruled by a set of three equations (continuity, Euler and 
Poisson) that in a ACDM model reads: 

^ + V-[(l + «5 m )v] = 0, 

P-+Hv +(vV)v = -V<£, 
or 

x/ 2 <t> =ln 2 n m s m , (i) 

where r is the conformal time. If we define the velocity di- 
vergence #(x, r) = V ■ v(x, t) and switch to Fourier space, 
the equations in (flj become: 



d5m(K t) 
8t 



+ 0(k,r) 



+ / d qd pfo(k - q - p)a(q, p)0(q, r)<5 m (p, r) = , 
9e{Kr) + ne(] i ,T) + ln 2 Q m (r)5 m (k,r) 



dr 



+ 



where 



J d\ d 3 p 5 D (k - q - p)/3(q, p)0(q, r)0(p, 



(2) 
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a(q,p) = 



(p + q) q 



/9(q,p) = 



(p + q) 2 pq 

2p 2 q 2 



(3) 



Equations ([2]) can be written in a compact form if we define a 
new time variable 77 = \n(D+ /D+i n ), being D+i„ the growth 
factor at an early epoch, and a doublet <p a (a = 1,2) 



V?i(k,7?) 



<M k ,f?) 

-0(k,r,)/{HU 



(4) 



with /+ = dlnD+/dlna. The velocity divergence is scaled 
such that it has the same dimension of the density contrast 
and in the linear regime —9(lc,rj)/(7if+) w <5 m (k, 77), i.e. 
<£i(k,0) = ¥>2(k, 0). The system is therefore 

d v v? a (k, 77) = — n a (, (77)^6 (k, 77) 

+e 7J 7 a6c (k, -p, -q)< / 9 b (p,77)( f o c (q,?7), (5) 

where sum over repeated indices and integration over 
repeated momenta are understood. The vertex function 
7a&c(k, p, q) (a, 6, c, = 1, 2) has only three non-vanishing el- 
ements 

1 



7121 (k, P, q) = - fo(k + p + q) a(p, q) , 
7222 (k, p, q) = (5 (k + p + q) /3(p, q) , 



(6) 



and 7112 (k, q, p) = 7121 (k, p, q), with 3d the Dirac-delta 
distribution. All the information about the cosmological 
model is contained in the matrix 



-1 



1 

2/s 



(7) 



that, in the following, will be considered as a constant ma- 



trix, approximating 



1. This ratio is indeed very close 



to unity for most of the history of the Universe. Making 
this approximation, one is modifying the behavior of the 
decaying mode, w hile the growin g one is left unaltered. It 
has been shown in lPietro ni (2008) that it affects the matter 
power spectrum at z — at a less than percent level up to 
k « 0.3/i Mpc" 1 . 

The power spectrum, defined by an ensemble average, 
in this notation is a 2 x 2 matrix 



(<A*(k; v)<fib(q.] V)) = M k + q)Pab(k; 77) ■ 
and the bispectrum, defined by 

(y>a(k; r/)</?&(q; T))<p c (p; 77)) 

= 8d(M + q + p)B abc (k, q, p; 77) , 



(8) 



(9) 



has 8 components. In the following, we will also consider a 
different-time two-point correlator, defined as 

{(p a Qt,r) a )<Pb(<l,r)b)} = S D (k + q)P ab (k; rj a ,rj b ) , (10) 
which obviously coincides with (|8} for 77a = 775 . 

2.2 The distribution of dark- matter haloes 

Let us consider a set of dark-matter haloes identified at a 
given redshift zid according to some predefined criterion. 
The material that forms the haloes can be traced back to its 
initial location in the linear overdensity field at z — > 00. 
We dub each of these regions as a proto-halo. In other 
words, a proto-halo is the Lagrangian region of space that 



will collapse to form a halo at redshift Zid. N-body simula- 
tions show that nearl y all proto-haloes are simply connected 
l|Porciani et al.l [2002T ) even though this property is not key 
to our analysis. 

Let us now follow the evolution of a proto-halo over cos- 
mic time in Eulerian space. Basically its shape and topol- 
ogy will be distorted (proto-haloes will first fragment into 
smaller substructures that will later merge to form the fi- 
nal halo) and its overall volume will be compressed while its 
centre of mass will move along a given trajectory determined 
by the mass density field via gravity. We focus our analysis 
onto this motion that connects the Lagrangian position of 
the proto-halo with the Eulerian location of the final halo. 

On scales much larger than the characteristic size of 
(and separation between) the proto-haloes, the density fluc- 
tuations traced by the centre-of-mass trajectories can be de- 
scribed with a continuous overdensity field (5h(x, r|zid). Note 
that while r labels conformal time along the trajectories, 
Zid is just a tag that identifies the halo population. Unlike 
real haloes that undergo merging, by construction proto- 
haloes always preserve their identity. Their total number is 
therefore conserved over time and we can write a continuity 
equation for their number density: 



<:)r 



+ V ■ [(1 + <5h)v h ] 







Here the proto-halo density and velocity fields should be in- 
tended as coarse grained on a scale of a few times the mean 
inter-halo separation (so as to suppress discreteness effects 
as proto-haloes are individually separate units). Strictly 
speaking, the smoothed velocity field does not obey the 
Euler-Poisson system in equation {TJ due to the presence 
of the non-linear term (v ■ V)v. In fact, the coarse grain- 
ing procedure introduces new terms in the fluid equations 
generated by the degrees of freedom one has integrated 
out, namely: a velocity dispersion term and a correction 
to the mean- field gravitational acceleration due to den- 
sity fluctuations on scales sma ller than the smoothing radius 
l|Buchert fc Dominguea l2005). On the other hand, it is rea- 
sonable to assume that the large-scale motion of the proto- 
haloes is generated by density fluctuations with wavelength 
larger than the characteristic halo size and is not influenced 
by perturbations with shorter wavelength. The very same 
assumption of neglecting the coupling to the small scales is 
routinely done when one writ es equation (flj) for the mass 
density field (see Section 3 in lBuchert fc Dominguezll2005l ) 
albeit adopting much narrower smoothing kernels than for 
the haloes. 

With the same spirit, in what follows, we will ignore the 
extra terms in the fluid equations generated by the coarse 
graining procedure. This is a working hypothesis which 
makes the problem mathematically treatable and whose ac- 
curacy will be tested by comparing our final results against 
high-resolution numerical simulations. We therefore write an 
Eulcr equation for the proto-halo fluid velocities 

<9v h 



+ Hv h + (v h • V)v h = -Vct> . 



(12) 



where the gravitational potential is the same as in equation 
(1). Note that if Vh matches v in the initial conditions then 
it will always do. On the contrary, any initial velocity bias 
between proto-halos and matter will be progressively erased 
by the gravitational acceleration. 



4 A. Elia, S. Kulkarni, C. Porciani, M. Pietroni, S. Matarrese 



Thus, given suitable initial conditions for 5^ and 
at r — > (i.e. a prescription for the Lagrangian biasing 
of proto-haloes) , we can in principle use equations and 
(|12[) to follow the clustering evolution of the proto-halo pop- 
ulation at all times. We are particularly interested in the so- 
lution of the fluid equations at the special time r that corre- 
sponds to Zid. In fact this solution has a particular physical 
meaning as it gives the density and velocity fields of the 
actual dark-matter haloes. 



2.3 Growth of matter and halo perturbations 

The system |TJ is now extended by the inclusion of eq. pip 
and eq. ((12}. We define a quadruplet ip a (a = 1, 2, 3, 4) 



/ *5i(k,?j) \ 
■MM) 



V 



e(Kv)/CHf+) 

Sh{k,v) 

9fc(k, !,)/(«/+) / 



(13) 



in such a way that eq. (JSJ still holds, but with indices running 
from 1 to 4. There are three more non-vanishing elements 
of the vertex 7343 (k, p, q) = 7334 (k, q, p) = 7121 (k, p, q) 
and 7444 (k, q, p) = 7222 (k, p, q), and the 4 x 4 ft matrix is 



ft = 



I 



\ 



\ 



(14) 



From the definitions © and @, with a = 1,2, 3, 4, we 
get a 4 x 4 matrix for power spectrum; in the following, 
we will focus on the matter power spectrum Pn and the 
matter-halo cross spectrum P13. 



3 INITIAL CONDITIONS 

In the previous section, we have presented a model that 
describes the non-linear evolution of the matter and halo 
density fields. Given suitable initial conditions, the formal 
equations we have introduced can be integrated numerically 
so that to compute the perturbative propagators and the 
TRG-evolved power spectra. The choice of the initial con- 
ditions therefore plays a very important role in our theory 
and will be the subject of this section. 



3.1 N-body simulation 

To gain insight into the properties of proto-haloes (and, 
later, to test our results at 2 = 0), we use one high- 
resolution N-body simulation extracted from the suite pre- 
sented by Pillepich et al. (2010). This contains 1024 3 dark- 
matter particles within a periodic cubic box with a side of 
Lbox = 1200ft -1 Mpc and follows the formation of structure 
in a ACDM model with Gaussian initial conditions and cos- 
mological parameters: h — 0.701, erg = 0.817, n s = 0.96, 
n m = 0.279, fit = 0.0462 and fi A = 0.721. 

We identify dark-matter haloes at 2 = using the 
friends-of-friends algorithm with a linking length equal to 
0.2 times the mean interparticle distance. We only consider 
haloes containing more than 100 particles (i.e. with mass 



M > 1.24 ■ lO 13 /!" 1 !'^) and we split them into four mass 
bins to keep track of their different clustering properties. The 
corresponding mass ranges and the total number of haloes 
in each bin are given in Table [1] along with an estimate of 
the highest wavevector up to which the fluid approximation 
for haloes holds. This value is determined by the number 
of haloes we require to be in a volume element to consider 
them as a fluid, and we set this number to 30. On smaller 
scales, our assumption breaks down, therefore we will look 
at results in the specified range, that, of course, decreases 
as the mass of the haloes increases. In the plots that will 
be shown in Section [5] the limit to which we can trust our 
model will be represented by vertical black dotted lines. 

Proto-haloes are identified by tracing the positions of 
the particles forming a halo at 2 = back to the linear 
density field. The centre of mass of each proto-halo is used as 
a proxy for its spatial location. Similarly, the mass weighted 
linear velocity gives the proto-halo velocity. 

Halo and proto-halo density and momentum fields are 
computed with the cloud-in-cell grid assignement using a 
512 3 mesh. Velocity fields are obtained by taking the ratio 
of the momentum and density distributions (preventively 
s moothed to pr e clude the existence of empty cells) as shown 
in IScoccimarrol (|2004T ). 

Power spectra have been computed using FFT. In order 
to avoid uncertain shot-noise corrections for the haloes, we 
only consider their cross spectra with the matter density 
field. 



3.2 Lagrangian halo bias 

Concerning the matter density, the initial conditions are 
given by linear theory which directly follows from the 
adopted cosmological model (transfer function) and the 
statistics of primordial perturbations (spectral index, Gaus- 
sianity). On the other hand, for the dark-matter halos, we 
can follow two different approaches: (i) extract the relevant 
information directly from the simulation or (ii) use a model 
for the Lagrangian bias of the halos. The latter option offers 
a number of advantages. First, it allows us to make gen- 
eral predictions independently of the simulation specifics. 
Second, it allows us to include halo bispectra in our for- 
malism (while it would be extremely demanding and time 
consuming to compute all possible triangular configurations 
from the simulation). For these reasons we will present be- 
low a model for the bias of the proto-haloes. Note, however, 
that any lack of accuracy of the adopted Lagrangian bias- 
ing scheme will propagate through the time evolution of our 
model and contribute to the imprecision of its final results. 
Therefore, in order to test the accuracy of our evolutionary 
equations alone, we will also extract initial conditions di- 
rectly from the simulation and compare the corresponding 
outcome of the evolution model with the statistics of the 
simulated haloes at 2 = 0. 

Let us consider the overdensity of proto-haloes in La- 
grangian space <5h(q) and the corresponding mass-density 
fluctuation <5(q). We assume that their Fourier transforms 
are linked by the expression: 



<54k) = (&i + &2-fc 2 )Mk), 



(15) 



which corresponds to a non -local relati o n in r eal space. This 
form was first proposed bv lMatsubaral l|l999l ) and describes 
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the clustering of linear density peaks |Desiacquesll2008l ). In 
this case, the initial conditions for P33 and P13 are: 



P 33 (k) = (&i + b 2 fc 2 ) 2 P 1 i(fc)e-' :2fl2 , 
Pa{k) = (b 1 + b 2 k 2 )P 11 (k)e- k2R2/2 . 



(16) 



where the exponential functions accounts for the finite size 
of the density peaks corresponding to a given halo mass. 
We find that the expression above accurately describes the 
cross-spectrum of proto-haloes and matter in our N-body 
simulation when bi, &2 and R are treated as fitting pa- 



rameters (see lElia et al.l 201 lj for fu r ther d etails), [j] This is 
not surprising as iLudlow fc Porcianil (|201fj| ) have shown that 
most proto-haloes include a density peak of the correspond- 
ing mass scale within their Lagrangian volume. In Table [1] 
we quote, for each halo-mass bin, the parameters 61, 62 and 
R that best fit the simulation data using eq. (|16[) where the 
linear matter pow er spectrum Pn i s computed using the 
CAMB online tool lLewis et alJ <|2000h . 

Whereas Gaussianity is a good approximation for the 
linear matter distribution, fluctuations in the halo counts are 
non-Gaussian even in the initial conditions. We can quantify 
their level of non-Gaussianity in terms of their auto and cross 
bispectra that, using equation (|15[) . can all be reduced to one 
of the following forms: 

P 3 33(ki,k 2 ,k 3 ) = (&i + &2fc 2 )(6i+b 2 fcl)(bi+6 2 fc 2 ) 

P m (ki,k 2 ,k 3 ) , 
P133 (ki, k 2 , k 3 ) = (61 + 6afca)(6i + 6 2 fef)P m (ki, k 2 , k 3 ) , 
S U 3(ki,ka,ka) = (&i + 6 2 fc 2 )B m (ki, k 2 , k 3 ) , (17) 

where the matter bispectrum P m (ki,k 2 ,k 3 ) is computed 
using the tree-level expression of the standard perturbation 
theory, 



'1 






k 2 \ I 2' 


2 + 


i 


\k 2 





Pm(ki,k 2 ,k 3 ) = 2 

Pn(fcl, «in)Pll(fe, z in ) + cycl. 

with Mi2 = ¥£■ 



(18) 



In order to solve our evolutionary equations, we need to 
know also the linear velocity field of proto-haloes. In princi- 
ple proto-haloes might not move with the same velocity as 
matter at the same location (at the very least they should 
match the mass velocity smoothed on the Lagrangian halo 
size). We model this effect assuming that proto-haloes are 
indeed related to linear density peaks which, as discussed 
above, gives a good descri ption of their cluster i ng pr oper- 
ties. In particular we follow iDesiacaues fc Shethl (|2010h who 
proposed a model for the peak velocities, which in Fourier 
space assumes the form 



h (k) = 1 



„2 
— J ft 



-k 2 R 2 /2 



m (k) = b v {k)6 m (k) , (19) 



with by the scale-dependent "velocity bias" and a n (n — 
0, 1) being the spectral moments of the matter power spec- 
trum defined as 



1 Note that, adopting a local-bias model, <5;,(q) = b± ■ <5 m (q) + 
62 ■ ^m(q)i provides a worse fit to simulation results. In this case, 
the model P13 lacks power for all but the smallest wavenumbers. 



1 

2^2 



dfcfc 2(n+1) Pii(fcK 



(20) 



In the equations above a Gaussian smoothing window of size 
R has been adopted to coarse grain the matter density and 
identify the peaks. Once the proto-halo mass is linked to 
R through M = (2tt) 3/2 p R 3 (with p the mean comovmg 
density of matter), this model has no free parameters, since 
the spectral moments are completely defined by Pn(fc). It 
follows that the initial conditions for P24 and P44 are 



P 44 (fc) 
P 24 (fc) 



&«(*) P 22 (fc) ; 
b v (k)P 22 (k), 



(21) 



and these expressions are in very good agreemen t with the 
spect ra computed from the N-body simulation (Eli a et al.l 
l201lf ). The corresponding values of o%ja\ are listed in Ta- 
ble [1] as a function of halo mass. Note that the velocity bias 
becomes more and more important on large scales with in- 
creasing the halo mass. 



4 LINEAR THEORY 

The lowest order approximation to the perturbation equa- 
tion ((51 consists in setting j aoc = 0. In this limit, the evolu- 
tion of the field from the initial time 77 = to a generic r\ is 
given by 

<p*{k;r,) = gab{rj)<Pb{k;G) , (22) 
where gab{v) i s the linear propagator, defined by the equation 



{Sabd v + flab) gbc{v) = SacSoiv) ) 



(23) 



with Sab the Kronecker delta. Solving eq. ([23} with causal 
boundary conditions (g a b(r i) = for 77 < 0, see e.g. 
ICrocce fc Scoccimarrd 120061 ) one gets 
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(24) 



with 6(n) Heavyside's step function. Notice that g a b(rf) — > 
8 a b as 77 —;• + . The first and second contributions repre- 
sent the standard growing an d decaying modes, respectively 
l|Crocce fc Scoccimarrd [20061 ) . The third and fourth contri- 
butions represent two new modes, decaying respectively as 
e -3?7/2 anc j e -v com p arec i with the growing one. To under- 
stand their physical effect we notice that an initial condition 
of the form 
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Bin 


Mass range 
(lO 13 M /ft) 


# haloes 


(Mpc" 1 h) 


6i 


b 2 

(Mpc 2 h~ 2 ) 


R 

(Mpc/i- 1 ) 


(Mpc 2 h~ 2 ) 


Bin 1 


1.24 - 1.8 


202948 


0.24 


7.28 ±0.38 


422 ± 102 


2.7 ±0.8 


9.1 


Bin 2 


1.8 - 3.4 


211305 


0.24 


14.2 ±0.4 


356 ± 84 


2.1 ± 0.7 


12.3 


Bin 3 


3.4 - 10 


150105 


0.22 


25.9 ±0.4 


708 ± 103 


2.9 ± 0.4 


20.5 


Bin 4 


> 10 


48985 


0.15 


66.2 ±1.3 


1025 ± 401 


3.5 ± 0.8 


50.9 



Table 1. Mass range and number of the haloes in the four bins. 



/ ¥>(k) \ 
V v«(k) / 

evolves (using eq. (|22|) ~) into </? a (k;?7) given by 



(25) 



tp ■ 



fv) + e~ v (-3tp + ifih + 2<Pv) 



(26) 



i.e. both the halo density and velocity fields relax to the 
matter ones as n — > oo (but at a different pace). Also note 
that in the absence of an initial density bias (i.e. 03 = (f>) but 
in the presence of an initial velocity bias (i.e. </>4 ^ </>), the 
linear dynamics quickly generates a transient density bias 
that vanishes at late times as e~ v — e -3 '"' 2 . The initial power 
spectrum at rj = 0, corresponding to the field configuration 
in eq. ((25}, is P® b (k), and it evolves forward in time as 

(27) 



PL,ab(k;r] a ,r] b ) = gac(j) a )g bd (j) b )P cd (k) . 



4.1 The importance of velocity bias 

It is interesting to assess the role of the velocity bias in the 
linear solution previously discussed. Assuming ip4 — ifi2 at all 
times, the linear propagator for the first three components 
ifii with i — 1, 2, 3 becomes 



9ab(rj) = 




and the third component of (|26[1 reduces to 
^3(k; J7 ) = ^(k) + e -' ) (^(k)-^(k)). 



(28) 



(29) 



This expresses the well known linear debiasmg between the 
halo and matter fields at late times, derived by iFrvl (|l996l ) 
for tracers that do not undergo merging and move solely 
under the influence of gravity. 

The corresponding halo-matter cross spectrum is 



L,13 



P ° 1+e -n { p0 3 _pO i)sPi 



(3) 



while keeping ip4 tfi2 one gets 



5(3) 
i,13 



(30) 



(31) 




Figure 1. A comparison between P14 (blue dashed line) and Pn 
(red solid line) in the four bins. The vertical black dotted lines 
represent the limit up to which we expect our fluid approxima- 
tion for haloes to work. A smoothing scale of R = 7Mpc//i has 
been used for P14. For a fair comparison with Pn, which is not 
smoothed, P14 has been redivided by the smoothing function. 



In Figure [T] we compare P{\ against P® 4 extracted from the 
N-body simulation, for the different mass bins. While the 
spectra agree well on very large scales (k < 0.05 h Mpc -1 ), 
they progressively depart for smaller scales. This is in line 
with the model introduced in Section 3. Note that the last 
term in eq. (|31[) vanishes in the initial conditions, reaches 
a minimum for 77 ~ 0.8 and it is suppressed at late times. 
We quantify its amplitude at z = in Figure O where we 
plot the ratio tl = Pl,X3/Pl 13 which ranges between — 3 
per cent, depending on halo mass and scale. This suggests 
that the effect of the velocity bias on the halo-matter cross 
spectrum is small for low redshifts. 



5 ANALYTICAL TREATMENT OF 
NON-LINEARITIES 

In this section we deal with the non-linear evolution of the 
matter-halo system. We first compute 1-loop corrections for 
the propagator and then perform the corresponding large-k 
resummation to all perturbative orders thus extending the 
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Figure 2. The ratio between Pl, 13 and P^ ( 3 : going from bottom 
to top, Bin 1 (red), Bin 2 (blue), Bin 3 (green) and Bin 4 (black). 



results presented by ICrocce fc Scoccimarrol (|2006t ) for the 
matter density field. Finally, we compute non-linear power 
spectra using the TRG approach. For simplicity we only con- 
sider the case with no velocity bias, which we have demon- 
strated to be accurate (at least in the linear regime) at low 
redshifts, where current observations are available. The 3x3 
ft matrix for this case is 



( 1 

3 
2 

V o 




(32) 



5.1 1-loop perturbation theory 

The 1-loop correction to the linear propagator 
ijCrocce fc ScoccimarrdfeoOot ) is given by 



dsi / ds2 I d qe 

u JO 



'l+"2 



g a c(n - si)7 cle (k, q, q - k)g ef (si - s 2 ) 
7/fed(k-q, q, -k)g d b(s2)PL,ih(q, si, s 2 ) , 

(33) 

(see appendix A for its explicit expressions). 

The 1-loop correction to the two-point correlator ()10[) 
is given by the sum of two contributions 



APi 6 (k; Va ,Vb) + AP'J (k;ri a ,r] b ), 



(34) 



which are also known in the literature as "Pis" and "P22", 
respectively 0. They are given by 

APai, (k;rj a ,Vb) = Ag ac (k;r/ a ) g bd (k-r] b ) P° d (k) 



2 They should not be confused with the P13 and P22 of our no- 
tation! 



AP ab (k;T] a ,Vb) 



with 



(35) 



g ac (k;r]a — s a ) gbd(k;r] b — sb) , (36) 



ds a \ ds b $cd(k; s a , s b ) 



$ cd (k;s a ,s b ) = 2e Sa+s " J d 3 q 

7 ce i(k, -q, q - k)P L>ef (q, s a ,s b ) 
Ph,ih{\k. - q|,s a , s b )"/ dfh (-\i, q, k - q) . 

(37) 

The explicit expressions for AP/jf are given in appendix A. 

5.2 Large-k resummation for the propagator 

In the large-fc limit, the 1-loop correction for the propaga- 
tor, eq. (|33[) . grows as fc 2 , and eventually dominates over 
the (fc-independent) linear propagator. Taking into account 
higher orders, the situation gets even worse. The 2-loop cor- 
rection grows as fc 4 , the 3-loop as fc 6 , and so on. This a 
manifestation of the perturbative expansion breakdown in 
cosmological PT, which appears not only in the computa- 
tion of the propagator, but also of the power spectrum, the 
bispectrum, and so on. However, for the case of the p ropa- 
gator, it was shown in ICrocce fc Scoccimarrd (|2006bl ) that 
the leading order corrections in the large fc and large 77 limit 
can be resummed at all orders in perturbation theory, giving 
a well-behaved propagator. 

The propagator G ac (k;?7) connects the initial correlators 
with the cross-correlations between final and initial field con- 
figurations, 



P ab (k; V ,0) = G ac (k; V )P° b (k) + APjI G (fc^,0) , 



(38) 



where the last term at the rhs comes from the initial non- 
Gaussianity of the matter and halo fields. At leading order, 
it is given by (see appendix [Cj 



APjT(M,0) 



dse s g a c(v ~ s)g d f(s)g eg (s) 



d qjcde(k, q,q k)B /s f,(q,k - q, -k) . 



(39) 



where B abc is the initial bispec trum at rj = (z — Zm) 
(see also iBernardeau et al. 2010). In Section loTTl we will use 
eq. (|38jl to assess the validity of different approximation 
schemes for the propagator. In the large-fc limit, G decays 



G ab (k; 77) = g a b{n) ex P - 



fcW" 



with 



,3 P\q) 
3 / dq — 



(40) 



(41) 



Therefore, at least in the case of the propagator, the bad 
ultraviolet behavior is just an artifact of the perturbative 
expansion, which, at any finite order, completely misses the 
nice -and physically motivated- G au ssian decay of eq . (I4TJ1) 
(see l Crocce fc Scoccimarr o 2006b and lMatarrese fc Pietronil 
I2007T for a detailed discussion). 
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Although the result (|40[) was obtained for the 2x2 
propagator of the matter density-velocity system, it holds, 
taking into account proper modifications, also when halos 
are included, i .e. for the 3x3 propagator con sidered in this 
Section. As in lCrocce fe Scoccimarrol (|2006bl ). we obtain an 
improved propagator interpolating between the 1-loop re- 
sult (eq. I33[) at low k and the Gaussian decay (|40[) (with a 
modified pre-factor) at high fc. The details of the derivation 
and the relevant formulae are given in appendix B. 



5.3 TRG 

Unlike the propagator, the power spectrum cannot be re- 
summed an alytically at large k. Different semi-ana l ytical 
procedures dCrocce fe Scoccimarrol |2006l . iMcDonaldl fepOTt 
iTaruva fc Hiramatsull2008r ) have been proposed to compute 
it in the mildly non-linear regime. In th is paper we will con- 
sider the TRG technique introduced in IPietronH (|200ct ) . 

Starting from eq. a hierarchy of differential equa- 
tions for the power spectrum, the bispectrum and higher 
order correlations is obtained. We choose to truncate it at 
the level of the trispectrum Q a bcd = 0, so that the equations 
for P a b and B a bc form a closed system 

d v P a b(k;rj) = 

- Q, ac (r))P cb {k; rj) - Q. bc {rf)Pac{k\ rj) 

+e n J d 3 q [7 a orf(k, -q, q - k) B bc d(k, -q, q - k; rj) 

+B acd (k, q, q - k; 77) J bcd (k, -q, q - k)] , 



d v B abc (k, q, q - k; 77) = 

-n ad (ri)B d b c (k, q, q - k; 77) 

-tt b d(r])Badc(k, q, q - k; 77) 

-n cd (rj)B abd (]<., -q, q - k; 77) 

+2e" [7„ de (k, q, q - k)P db (q; 77)P ec (|k - q|; 77) 

+7bde(-q, q - k, k)P dc (|k - q|; i])P ea (k- 77) 

+ 7cde (q - k, k, -q)P da {k; T])P eb {q; 77)] , (42) 

which integrated gives the power spectra at any redshift 
and for any momentum scale. The system (|42p consists of 
coupled differential equations which are solved numerically, 
starting from given initial conditions, i.e. P a b(k; rji) and 
B abc {k, q, q k; 77,). 

From eq. (|32jl . we can observe that f2i3 and Q23 are zero, 
which means that the evolution of 8 m and 9 is not modi- 
fied, with respect to the original TRG formulation, by the 
presence of Sh, as it is expected. 



6 RESULTS 
6.1 Propagator 

In order to assess the validity of our analytical approach, we 
compare our results for the resummed propagator against 
the simulation; to this end, we consider the relation in 
eq. (|38[) . In particular we choose the indices a — 3, b — 
1, so that we can check the components related to the 
haloes that were not present in the original formalism by 
ICrocce fc Scoccimarrol l|2006bf ) and, at the same time, we 



0.01 = 



0.001 



Bin 1; 




Bin 2: 



0.01 
k (h Mpc" 1 ) 

Figure 3. Cross spectrum between the halo density at z = and 
the matter density at z = 50. The outcome of the N-body simu- 
lation (black points with error bars) is compared against linear- 
theory PpL (blue dashed line) and resummed result Ppji (green 
solid line) . The red dotted line shows the effect of neglecting the 
non-Gaussian term, i.e. Ppjj — AP^ 5 . 

The vertical black dotted lines represent the limit up to which we 
expect our fluid approximation for the haloes to work. 

do not have to deal with the shot-noise problem. We ex- 
tract the cross spectra from the simulation and compare 
them against those obtained both using linear theory prop- 
agators P PL = QziPii + S32-P21 + ff3i-P3°i and resummed 
propagators P PR = G 3 iP?i + G32P21 + Gsiffi + AP 3 T°; 
the result is shown in Figure [3] We note that the linear 
model severely overpredicts the two-time cross spectrum for 
k > 0.05 h Mpc -1 . It is evident that the resummed the- 
ory improves considerably upon the linear one, and agrees 
with the simulation within 10 per cent accuracy up to the 
scale where the fluid approximation holds. We include in 
AP£{ G the effect of the initial non-Gaussianity of the halo 
field via its initial bispectra, computed as in eq. (|17|l . It 
turns out that the components giving a non-vanishing ef- 
fect are of the -B113 type. Their contribution is suppressed 
by a D+(zin)/ D+(z = 0) factor with respect to that of an 
hypothetical primordial non-Gaussianity in the matter field 
(which we do not consider here). Therefore, the effect is of 
modest entity but, nevertheless, it improves the agreement 
with the simulation with respect to the case in which it is 
neglected. 



6.2 Power spectrum 

The TRG equations presented above are integrated numeri- 
cally starting from the initial conditions discussed in section 
[3] As a first step, we set all the initial bispectra to zero In 
Figure [4] we show a comparison between the halo-matter 
cross spectra extracted from the N-body simulation and 
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Figure 4. The cross spectrum between matter and halo distri- 
bution at z = is shown in the four bins; the black dots with 
error bars represent the simulation, the blue dashed line is linear 
theory, the violet dot-dashed line is 1-loop and the red solid line 
is TRG. The vertical black dotted lines represent the limit up to 
which we expect our fluid approximation for the haloes to work. 

the results of TRG, one-loop and linear theory. In the first 
three bins, corresponding to lower halo masses, linear the- 
ory overpredicts the power on mildly non-linear scales; note 
that this departure arises on larger scales compared to the 
matter auto spectrum (not shown in the figure). The over- 
prediction of linear theory is cured by the one-loop power 
spectrum only on very large scale, while the TRG manages 
to correct it up to a smaller scale, before starting to fail. 
The fourth bin, though, displays a different behaviour: lin- 
ear theory lacks of power on small scales, and neither the 
1-loop correction nor the TRG are much more accurate. This 
might originate from the fact that very massive haloes are 
large and rare in the initial conditions, therefore less suited 
for the fluid approximation. Moreover, they also display the 
strongest velocity bias, which we are neglecting in our cur- 
rent non-linear treatment. A more quantitative analysis is 
presented in Figure [SI where we plot the ratios between the 
spectra from the simulation and the theoretical results. The 
TRG gives a cross spectrum within 5% accuracy at least up 
to k — 0.1 ft Mpc -1 (barring bin 4), while linear theory does 
so up to k = 0.05 ft Mpc -1 . 

As previously pointed out, however, the halo-density field is 
not Gaussian at Zi„ = 50. We use expressions l|17p to ac- 
count for initial bispectra in the TRG method (they are not 
present in the linear theory and at 1-loop level) . Figure [5j 
illustrates that this indeed produces a slight improvement 
in the agreement with the simulation. The correction re- 
sulting from the introduction of the initial bispectra turns 
out to be quite small. The reason is that their contribu- 
tion to the final cross spectrum P13 is suppressed. To un- 
derstand why, we can use perturbation theory; first, let us 



investigate the case of Ph. In the one- loop computation, 
the initial Pu contributes to the final one with a "weight" 
of (D+(z — 0)/ D+(zin)) 2 = e 2v , as it also happens in lin- 
ear theory. The initial matter bispectrum Bin, instead,has 
a weight of e v , so it is suppressed by a factor of e~ v q If 
we now consider the haloes, we have showed in (|28p that 
there is a new decaying mode, responsible for the linear de- 
biasing. This new mode, that involves only the halo field, 
carries an extra e~ v suppression factor. We can now rank 
the contributions to P13 according to their relevance: 

0) P11 

(ii) P13, Bin 

(iii) P33, B113 

(iv) B133 

(v) B333. 

Each item is suppressed by a factor of e~ v with respect to 
the previous one. We can see that only Bin has some rele- 
vance, while the other terms are highly suppressed. Even if 
the reasoning was based on perturbation theory, it is valid 
also for TRG, at a qualitative level. Incidentally, the fact 
that Pn is the most relevant contribution is another evi- 
dence of debiasing. 

We can now address the effect of the truncation we perform 
in the TRG, namely considering the trispectrum Q = 0. 
First of all, the matter trispectrum Qim can be neglected 
in the range of scales under consideration, as one can con- 
cl ude from the c omparison between TRG and simulations 
in lPietronil (|2008l ). The contribution from initial mixed (i.e. 
matter-halo) or pure halo trispectra is furtherly suppressed 
with respect to Qim by extra e~ v factors, for the same 
reason as above. However, the trispectrum has its own time 
evolution as well, and one might argue that it becomes more 
relevant for z < Zi n ; this seems to be not the case, because 
enlarging the TRG truncation scheme by including the run- 
ning of the trispectrum gives a contribution to the power 
spectrum which is at least of two-loop order and it is cer- 
tainly subdominant in the scales we are considering. 
From Section 5 onward, we neglected any velocity bias, since 
this approximation proved to be accurate enough at z — 
(in linear theory). As a further check, it is interesting to ob- 
serve from Figure that the TRG is able to give a better 
prediction than full linear theory in eq. (|31[) , even neglecting 
the velocity bias that the linear theory accounts for. 
We can now look at the model predictions for the halo bias, 
defined as the ratio P13/P11. This quantity is plotted as a 
function of wavenumber in Figure [7] While the linear-theory 
bias always increases with scale, irrespectively of halo mass, 
the TRG result closely follows the scale dependence of b(k) 
seen in the simulation for the first two bins. It also gives 
a nearly constant bias for the third bin, as the simulation 
does, even though with a slightly lower value. However, the 
linear model performs better in the last bin where the bias 
in the simulation increases with k. 

It is also interesting to look at the results displayed 
in a different way; we can investigate the evolution of the 
cross spectrum from the initial conditions to today and the 



3 For the most experienced readers, this happens because one of 
the two vertices carrying the e* 7 factor is replaced by the bispec- 
trum. 
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Figure 5. The ratio between the spectrum from simulation and: 
linear theory (blue squares), TRG without bispectra (red trian- 
gles), TRG with bispectra (green circles). The vertical black dot- 
ted lines represent the limit up to which we expect our fluid ap- 
proximation for haloes to work. The shaded area marks the 5 per 
cent accuracy interval. 



Figure 7. The effective bias b = Pmh/Pm as a function of the 
wavevector in the four bins: the black circles with error bars are 
from the simulation, the blue dashed line represents linear theory 
and the red solid line TRG. The vertical black dotted lines repre- 
sent the limit up to which we expect our fluid approximation for 
the haloes to work. 
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Figure 6. The ratio between the cross spectrum from simulation 
and TRG with bispectra (green) and linear theory with the in- 
clusion of velocity bias (magenta). The vertical black dotted lines 
represent the limit up to which we expect our fluid approximation 
for haloes to work. 



evolution of the bias as well. In Figures [8] and [9] we plot, 
respectively, 



P mh (z = 0) 
P mh (z = 50) 



and 



b(z = 0) 
b{z = 50) 



(43) 



Again, our model is able to match accurately the trend of 
the simulation and to improve upon linear theory, excluding 
bin 4. A key feature of the linear solution in eq. ()26[) is the 
debiasing between halo and matter distributions with time. 
It is worth noting from Figure [9] that this effect is stronger 
for high-mass haloes, but constant on all the scales, while 
for low-mass haloes it is weaker on large scales and presents 
a strong fc-dependence. 



7 CONCLUSIONS 

We have presented a novel approach to modeling the clus- 
tering of dark-matter haloes on mildly non-linear scales. 
This follows the motion of the regions that will collapse to 
form haloes (that we dub proto-haloes) . Since the number 
of proto-haloes is conserved over time, for sufficiently large 
scales (k < 0.2 h Mpc" 1 ), we can write a set of fluid equa- 
tions that govern their evolution under the effect of gravity, 
which couples perturbations in the halo and matter den- 
sity fields. We provide analytical solutions for the linearized 
equations and 1-loop perturbative corrections for the halo 
and matter power spectra. For the propagator, quantifying 
the memory of the density and velocity fields to their ini- 
tial conditions, we also perform a resummation of perturba- 
tive corrections. Finally, for the power spectrum we compute 
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Figure 8. The ratio between the halo-matter cross spectrum at 
z = and z = 50 for the four mass bins: symbols represent the 
simulation, dashed lines the linear theory and solid lines TRG. 
From top to bottom, we have Bin 1 in red, Bin2 in blue, Bin 3 in 
green and Bin 4 in black. 
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the non-linear evolution using a semi-analytical procedure, 
namely an extension of the time renormalization group. 

The initial conditions for our evolutionary equations are 
specified adopting a Lagrangian bias model, originally devel- 
oped to describe the clustering and motion of linear density 
peaks. We fix the parameters of the model so that to re- 
produce the distribution of proto-haloes in a high-resolution 
N-body simulation at z = 50. We use the same simulation 
to test the predictions of our model at z — 0. 

Our main results can be summarized as follows: 

• Independently of the initial conditions, in the linear so- 
lution the halo density and velocity fields asymptotically 
match the corresponding matter fields at late times. This 
'debiasing' develops at a different rate for the density and 
the velocity, being faster for the latter. 

• Even if there is no initial density bias, the presence of a 
velocity bias generates a transient density bias that vanishes 
at late times. 

• Neglecting any initial velocity bias alters the linear pre- 
dictions for the halo- matter cross spectrum at redshift z — 
only by less than 3 per cent, for k < 0.3/iMpc~ . This pro- 
vides us with the motivation to ignore the velocity bias in 
the non-linear analysis. 

• Unlike its linear counterpart, the resummed propagator 
is in good agreement with the N-body simulation, indepen- 
dently of halo mass. 

• The halo-matter cross spectrum predicted by the TRG 
is accurate to 5 per cent up to k ~ 0.1 h Mpc -1 for a broad 
range of halo masses. This does not hold for very massive 
haloes (M > 10 14 /i _1 Mq), that have low number density 
and high initial velocity bias, for which discreteness effects 
are more important. 

• The TRG result improves upon both linear theory and 
1-loop corrections. Its performance is slightly enhanced ac- 
counting for the initial non-Gaussianity of the halo distribu- 
tion. 

• For low halo masses our model accurately describes the 
scale-dependent bias measured in the simulation at z — 0. 
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Figure 9. The ratio between the bias at z = and z = 50 for the 
four mass bins: symbols represent the simulation, dashed lines the 
linear theory and solid lines TRG. From top to bottom, we have 
Bin 1 in red, Bin2 in blue, Bin 3 in green and Bin 4 in black. 
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APPENDIX A: COMPLETE ANALYTIC EXPRESSIONS OF THE 1-LOOP PROPAGATOR AND 
POWER SPECTRUM 

Linear power spectra: 

P L ,ii(k;r),ri) = Pn(fc) = P L , 12 (k; 11, r/) = P L ,22(k;i],i]) , 
P L , 13 (k;V,v) = e-"((e , '-l)i^(fc) + f&(fc)), 

P £ , 33 (M,r?) = c- 2 "(( e ' ) -l) 2 P 1 1 (fc) + P3 3(fc) + 2(e' ) -l)P 1 3(fc)). (Al) 



Al Propagator 

First we report the integrands of the one-loop corrections of the components we are interested in. The linear power spectrum 
can be split as 

Ph.ci.cj (qv, Si, Sj) = P°(qi)Uc t , Cj + AP CitCj (qi;Si,Sj) , (A2) 
where 

U CilCj = 1 for any a, Cj , (A3) 
while 

AP Ci , Cj (ft; Si, Sj) / only if a or Cj ; = 3 . (A4) 

Looking at eq. (|33|) we can consistently split the one-loop correction to the propagator as Ag a b = A^g a b + A^g a b and we 
denote them as 5^ % 'g a h, so that A^'g a h — J dqd^'g a h, for i = 1, 2. 



«5S 



li 



(e" - l) 2 Pn (q) [4 (6fc 7 g - 79fc 5 g 3 + 50fc 3 g 5 - 21fcg 7 ) - 3 (fc 2 - g 2 ) 3 (2fc 2 + 7g 2 ) log 



840fc 3 g 3 

5 (1) S3i = <5si 



7ii ■ 



<5 (1) ff32 = |5ffll , 

tfWflas = -| e -"( e '-l)V7rPf x («) : 



.931 



3" (^i° 3 (?) - P?i(q)) (4 (9fc 5 g + 24fc 3 g 3 - 9fcg 5 ) - 18 (fc 2 - g 2 ) 3 log (|f±||)) 



280fc 3 g 



jc(2) 2 ( 2 ) 

5 (2) ff33 = (A5) 
A2 1-loop PS 

From eq. ()37|) . one can see that the one-loop corrections require integrals over ds\, ds2 and d 3 g = g 2 sin 6 dg dd d<j>\ the 
integration over q and over /1 = cos cannot be performed analytically. The following expressions are therefore integral 
kernels, denoted by S qlJ ,P (the factor of g 2 is already taken into account). If cos 9 = /j, = -jf 3 and \k — q\ — \JkP- — 2k[iq + g 2 , 

TrfcVPfUg) [7fc M + (3 - 10 M 2 ) g] 2 P° (|fc - g|) 



A P 77 



49 (A; 2 - 2fc/jg + g 



21-' 



(A6) 



u q^ r 13 



= ^^ol^t/ 3 N fc - (|fc - 1\) + i£(9)7/i - 2fc M g + g 2 ) P° 3 (|fc - ,|) 

49 («r — 2fc/ig + q z ) 

+kP? 1 (q)P? 1 (\k-q\) (7kn (e" - 1) - g ((10 M 2 - 3) c" - 14 M 2 + 7))] . (A7) 

= .., |2 ; 2 , 2 {98 W (fc - /ig) (fc 2 - 2fc M g + g 2 ) P? 3 (q)P? 3 (\k - q\) 

49 (k z — 2/c/ig + q z ) 

+7q(k - Qfc - g|) [2fcP° 3 (<z) [7fc M (e" - 1) - 4 ((10 M 2 - 3) e" - 14 M 2 + 7)] + 7qP^(q)(k - M )] 

+fc 2 P 1 ° 1 (g)P 1 ° a (|fc - g|) [g ((10 M 2 - 3) e" - 14 M 2 + 7) - 7fc At (e" - 1)] 2 
+49 M 2 P 1 ° 1 (g) (fc 2 - 2k m + g 2 ) 2 P 3 ° 3 (|fc - q\) 

+14fc A tP 1 ° 1 (g) (fc 2 - 2km + q 2 ) [7/c/J (e" - 1) - g ((lO^ 2 - 3) e" - 14 M 2 + 7)] P? 3 (|fc - g|)} . (A8) 
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APPENDIX B: THE RESUMMED PROPAGATOR 

As it was discussed in ICrocce &z Scoccimarr 3 l|2006bl ). the leading contribution in fc 2 exp(2n) to the propagator at a generic 
n-loop order contains a chain of propagators and vertices of the form 

gab ± {V - Sl)7i>l ci ai (k, -qi, Ql - k)# ai f, 2 (si - S 2 ) • • • Sa 2n _i 6 2 „(s2n-l - S 2 n)7i, 2 „ c 2 „ a 2 „ (k + q 2 „ , ~q 2 n , -k)g Q2n b (s 2n ) . (Bl) 

The d indices have to be contracted in all possible pairings, by inserting n linear power spectra, each of the form 

fo(qi + <te)Pi,,c i ,c j (qi;si,Sj) . (B2) 

The n-loop contribution is obtained by multiplying by exp(^ 2 ™ 1 Si) and by the appropriate combinatoric factor, and then by 
integrating over 

n£i (T ds z j d 3 g?j . (B3) 

Recall eqs. (|A2HA4[) : if all the insertions are of the P°(qi)U CilC type then the resummation goes exactly as in the standard 
case. Indeed 

,l3_ k 

2 q 2 

in the fc 3> q limit. Besides the explicit form for the vertices (see Section 1), we have used the composition property of the 
propagators 

g a b(v - si)gbc(si - s 2 ) = g ac (ri - s 2 ) , (B5) 

and have defined u = (1, 1, 1). Since each vertex, contracted byati c vector becomes proportional to a delta function in its 
first and third index, the chain of propagators composes up to a single propagator g a b(r\) and the time integral can be easily 
performed. The momentum integrals factorize into n integrals of the type 

dy^M^^tV 2 , (B6) 

T 

where a 2 has been defined in (141[) . Using the appropriate combinatoric factors, the leading contribution to the propagator at 
n-loop, in the large momentum limit, when all the power spectrum insertions are of the P° (qi)U Ci ,cj type is 



E/ci,a f flo6 1 (»7-sx)7i>ie 1 o 1 (k,-qi ) qi -k)g ai b 2 (si - s 2 ) ->■ u C] g ab2 (n - s 2 )-— 5- , (B4) 



1 



-k 2 a 2 ^l 



g* b (v) , (B7) 



which resums to g a b(fl) exp[— fc 2 a 2 ^ e ^ ]. 

As for the insertions including the AP Ci]Cj (qi\ Si, Sj) contribution to the linear power spectrum, the important point to 
realize is that, due to (|A4|) and to the structure of the linear propagator and the vertices, a chain like (|B1[) can be contracted by 
at most one AP CijCj . (qij Si, Sj) power spectrum, the remaining ones being of the type P (qi)U Ci , Cj ■ Moreover, these insertions 
only contribute to G31 and Gs2 - Therefore, at n-th order, we have only two types of contributions: those with all P° (qi)U Ci , Cj 
insertions, giving (IB7|I . and those with one AP CijC . insertions and n— 1 P°(qi)U Ci ,c- ones, which can also be resummed in the 
large k limit. 

As a consequence, the complete resummed propagator in the large momentum limit is 

G ab (k; V ) = (g ab ( V ) + 5 a3 f b A (2) g 3 i (fc; V)) exp[-fc 2 a 2 (e " ~ 1} ' ] , (B8) 
with f b = (1, 2/3, 0). 

In order to have an expression for the propagator valid at any k, one can proceed as in IjCrocce fc Scoccimarr 3 (|2006bl ^ 
and interpolate between the large k limit above and the 1-loop result at low k. This can be done, for instance for Gn, starting 
from its 1-loop expression 



(B9) 



511(77) + Agu(k;ri) ~ 311(77) + ^A gil (k;r))j 
where the above approximation is exact in the large n limit. Since 

|Affn(fc;n) ^ -fcV (£ " ~ 1)2 , for large fc, (BIO) 

the required interpolation is given by 
"5 



Gn(fc;n) =pn (»?) exp 



3 A 3 ii(fc;n) 



(Bll) 



Proceeding analogously for the other components, we have: 
Gi 2 (fc;«) = #i 2 (r?)exp 



5 

-Agii(fc;n) 
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G 2 2(fc; 77) 

Gai(M) 
(?3i(fc; 77) 
G 3 2(k;ri) 



-Ag 22 (k;rj) 



g 2 2(v)exp 
32i(?7)exp 
(33i(»7) + A (2) g 3 i(fc;?7)) exp 



-Ag 2 2(k;ri) 



-Agu(k;r)) 



332(7;) + -A (2) ff 3 i(fc;77) ] c xp 



r Agii(fe;77) 



G33(fc;77) = 


333 (77) exp [e^A^^ 


where we have used the identities 


Agi 2 (k;r]) 


2 

= gAffu(fc;r7), 


Ag 2 i(k;i]) 


3 

= -Ag 2 2(k;ri) , 


A (1) 33i(fe;7?) 


= Agn(k;ri) 


A (1) 33 2 (fc;7;) 


= |A (1) ff3 i(fc;r?) 


A (2) 33 2 (fc;7;) 


= |A«ff3i(fc;r?). 



(B12) 



(B13) 



APPENDIX C: TAKING INTO ACCOUNT THE INITIAL HALO NON-GAUSSIANITY. 

While at Zi n the matter field can be considered gaussian with high accuracy, the same not necessarily holds for haloes. An 
ini tial non-Gaus s ianity for the equal-time power spectrum can be easily incorporated in the TRG formalism as discussed 
m iBartolo et all (|20ld ). and has been done in sect. 6.2. In order to obtain the effect of an initial non-Gaussianity on the 
cross-correlator at different times considered in eq. (|38[1 . it suffices to compute the 0(7) correction to the linear evolution of 
the field, by inserting eq. (|22p in (0, to get 



<Pa ) ( k > r i) = J dsg ac (rj - s)e J d 3 q j cde (k, -q,q-k)gdf(s)(pf(q;0)g eg (s)ip g (k -q;0). 
The cross-correlator 

(^(Mtotk'jO)), 

then includes the non-Gaussian expectation value 

(?/(q; 0)<p a (k - q; 0)^(k'; 0)) = S D (k + k')B/ 9 fc(q, k - q, -k) , 

and therefore gives eq. (J39J) . 
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